The lifetime of Bc-meson and some relevant problems 



Chao-Hsi Chang^'^, Shao-Long Chen'^, Tai-Fu Feng^ and Xue-Qian Li^'^'^ 
1 CCAST (World Laboratory), P.O. Box 8730, Beijing 100080, China 
Institute of Theoretical Physics, Academia Sinica, P.O. Box 2735, Beijing 100080, China 
^ Department of Physics, Nankai University, Tianjin 300071, China 



Abstract 

The lifetime of the Bc-uieson is estimated with consistent considerations 
on all of the heavy mesons {B^, B^, Bg, D^, D^Ds) and the double heavy 
meson Be- In the estimate, the framework, where the non-spectator ef- 
fects for nonleptonic decays are taken into account properly, is adopted, 
and the parameters needed to be fixed are treated carefully and deter- 
mined by fitting the available data. The bound-state effects in it are also 
considered. We find that in decays of the meson B^ the QCD correc- 
tion terms of the penguin diagrams and the main component terms ciOi, 
C2O2 of the effective interaction Lagrangian have direct interference that 
causes an enhancement about 3 ~ 4% in the total width of the Be meson. 

14.40.Nd, 14.40. Lb, 12.39.Hg, 12.38.Lg 
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I. INTRODUCTION 



Recently the meson Be has been observed by CDF collaboration at Tevatron [|^], that 
careful studies of the meson Be are motivated with the fresh reason. 

It is known that B^ is a meson of the ground state of a double heavy flavored quark- 
antiquark system (there is no light flavored quark as a valence being involved), and its 
decay must be either through the decays of each component (a heavy flavor), or through 
the annihilation of the two components (two heavy flavors), that is very different from B 
and D mesons. Of the decays, the contribution from each component (6-quark or c-quark) 
to the total width happens to be comparable each other so in the future experiments it is 
accessible that to have precise measurements on each component. Namely with the meson 
of Be alone, one may investigate the two different heavy flavors simultaneously. Especially, 
certain decay mechanisms play a similar role in D-decays, -B-decays and -Bc-decays so some 
parameters appearing in Bc-decays should be the same as those in B-decays or in D-decays, 
therefore, when one estimates the -Bc-decay one may use the experimental available data of 
/^-decays and -B-decays as input phenomenologically to determine them under the consistent 
considerations. Obviously in this way the estimate for B^ lifetime should be comparatively 
reliable. 

The meson Be certainly is an independent complement to S-mesons and D-mesons for 
studying the two heavy flavor b and c decays. Furthermore, if one carries on a comparitive 
study of the two heavy flavors, it has unique advantages. In this sense, i?c-nieson will offer 
an extra interesting and unique laboratory for the heavy flavor decay studies. 

There have been quite a lot of studies on the lifetime of D and B mesons and the meson 
Be as well [0-0 The reason in part is that for the lifetime it is comparatively 'easy'. Due 
to the duality for quarks and hadrons: 

h^9j){(li,9j\ = \hk){hk\ 

ij k 
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where hk, Qi and gj denote hadrons, quarks and gluons respectively, the optical theorem may 
be set on the level of hadrons or the level of quark-gluons, an inclusive processes of hadrons 
can be turned onto a quark level instead. In general, the problem for evaluating a decay 
rate of a hadron is hard, because the relevant hadron matrix element cannot be handled 
reliably. The matrix element contains non-perturbative QCD effects, so one cannot compute 
them very satisfactorily based on an existent underlying-theoryQ. The optical theorem can 
be applied for evaluating lifetimes and certain inclusive processes, so the problem can be 
'solved' in part: the non-perturbative effects are summed by the theorem, generally only in 
the initial state are still needed be handled. Thus in the studies of the lifetime and/or the 
inclusive process as well, one may focus the efforts mainly on the decay mechanisms. 

For the estimate of the lifetimes and inclusive decays, the effective Lagrangian for weak 
decays with QCD corrections should be known 0-^. With the effective Lagrangian for c 
and b decays, phenomenological analyses of D meson lifetimes and S-meson lifetimes have 
been made For all the heavy meson decays the contributions can be decomposed into 
three categories: the dominant one i.e. the direct decay of the heavy quark while the light 
quark remains as a spectator (this contribution is very sensitive to the heavy quark mass 
i.e. proportional to Mq); the non-spectator one from W-annihilation (or exchange) (WA or 
WE); and the one from the Pauli interference (PI) [H. The parameters which are needed to 
be fixed are the quark masses, the matrix element {0\Ji^i5\M(^Bc,b,d)) relating to the decay 
constant, and the relevant non-factorizability parameters etc as well [^. The range of all 
the parameters are known, but their precise values are not well calculable. The better way 
is to fix them phenomenologically by fitting data. 

To have a better estimate of the lifetime of Be than before, we will take a 'consistent' 

-^In principle the lattice gauge simulation may deal with the non-perturbative effects as well as 
one wishes, but in practice the computer ability now still is at quite sizable 'distance' to obtain 
sufficiently accurate results for calculating such hadron matrix elements. 
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view of the parameters appearing in estimating the hfetimes for all of the heavy mesons 
D, Dg, B, Bg and those in estimating that for the meson B^. Namely to estimate the lifetime 
of the meson Be with the parameters fixed by phenomenologically fitting the available data 
for the other heavy mesons. We also try to discuss some uncertainties of the estimate in the 
paper. 

In literatures, the charm quark mass rric appearing in the estimate for D and B decays 



takes different values [pip!5|. We think it is reasonable, because the mass appears in different 
situations: in the initial state for D-decays but in the final state for B-decays. In general, 
for the quark (antiquark) in the parent meson of a concerned decay mode, its mass should 
take its 'pole' value if the bound-state effects are ignored, whereas, for the masses of the 
product quarks (antiquarks) in the final state of an inclusive process, it is more reasonable 
to take relevant running masses and the running energy-scale should be the mass of the 
decaying quark (or mesons for WA and PI). Anyhow, this problem is somewhat subtle. In 
the earlier estimations for inclusive processes, the quark decays are considered only as if 
the quark is 'free'. However, some authors have pointed out that the bound-state effects 
on the effective mass of the heavy quark should be taken into account. Namely the heavy 
quark effective mass, appearing in the formulation, should deviate from the pole value by an 



amount to correspond to the binding energy IXGHlgl . In our work, we also pay attention on 



the effects and use a parametrization which is a bit different from that of [|T^ to account for 



the bound-state effects on the mass of decaying quark (see the context below for details). 



The relation between the pole mass and the running MS mass is used many times in 
our estimation, so we present it up to one loop level here precisely for convenience. It reads 

m = m[m){l + - j, (Ij 

O TC 

where the running coupling constant at 1-loop level is 



127r 

(33 - 2nf)ln^ 
with as{m1) = 0.118 fl^. The running mass runs as 



= t:::^ — :r^r^^ (2) 



For the lifetimes of B and D mesons, the contributions from penguin terms of the effective 
Lagrangian generally are not important because of smallness of their coefficients C3 ~ cg. 
But as pointed out in |]19| , the penguin contributions to the charmless decays of B- mesons 
are not negligible. The reason is that for those modes the main contributions (since this 
parts are not zero, if we return back to the tree level, thus we will call them as 'tree parts' 
as in the most literature) suffer a cancelation (ci + C2/Nc) or (c2 + Ci/Nc), and the 'tree 
part' CiOi + C2O2 does not contribute in addition, thus the penguin contributions become 
important. 

As for Be meson, the problem has not been investigated very carefully. For instance, 
in the earlier paper 0, the lifetime of Be was estimated where the bound state effect was 
carefully handled in terms of the Bethe-Salpeter equation, but in the effective Lagrangian 



the penguin contributions were ignored. Recently, Beneke and Buchalla |Tl[ also presented 
an evaluation of the B^ lifetime, where they also ignored the penguins'. For the spectator 
mechanism, the contribution from the penguin terms in B decays has been estimated by 
Bagan et al. 0, and their results show that at most it can give rise to a few thousandths 
of changes, so in general it can be neglected. However, for the WA and PI terms, the 
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operators induced by the penguin diagrams are CiOi which contain terms {sibi){cjCj) and 

i=3 

{sibj){cjCi), where i,j are color indices, so in some Be decays they can interfere with the 
'tree part' 

^efT^ = VcbVe^lciCLlf.biSLYcL + C2S L'^^biClY Cl]- 

Therefore, in B— and decays, the contributions relating to penguins can be proportional 
only to or \c*Cj\{i,j = 3 ~ 6), so they are small whereas, in some Be decays there 
may exist \c*Ci\ and \c*C2\ terms, which are not so small, and such interference may bring 
up a few percents of corrections in its lifetime. Indeed our numerical results show the fact 
that the interference can make a change in lifetime of Be so large as 3 ~ 4% of the total. 



In the paper, including the penguin contributions and fixing the parameters with a 
'consistent' view to fit the data of B, Bg and Dg decays, we re-estimate the hfetimes, the 
branching ratios of the semileptonic and pure leptonic decays of the meson Be- We expect 
to gain more information on Be about QCD and decay mechanisms. 

The paper is organized as follows: after the introduction, we present the formulation in 
Sec. II, and give the numerical results and the concerned phenomenological parameters in 
Sec.III, then put conclusions and discussions in the last section. For convenience, we collect 
some useful formulars in Appendix A. 

II. FORMULATION 

In this section we will describe the different mechanisms to the lifetimes for the mesons 
B and B^ etc., and present the formulas for later numerical calculations. 

A. The spectator components and the contributions from 6 or c decays 

With quark-hadron duality and the optical theorem, the 'full' inclusive decay width 
(the lifetime) of a heavy hadron Hq (containing a heavy quark Q = b,c) is related to the 
absorptive part of the forward matrix element of the transition operator T. 

r{HQ ^X)^ -^Im f d'x{HQ\f\HQ) = -^{HqIFIHq), (4) 

where 

f' = TO£e//(x),£e//(0)}, 

and Ceff is the relevant effective weak Lagrangian which is responsible for the decay. For the 
concerned final state X with designated quark- antiquark combination, up-to order l/mg we 
have: 



6 



- A-) = ^IV(CKM)P{4{H,\QQ\H„) + ,x (£o!*:l_GQ|£«) 

+ E-J.™^#^=™ + 0(1A4)}. (5) 

Here only the heavy quark (6, c quark) decays are concerned. In the spectator components 
of the decays, for the heavy meson decays, the hght flavor in the heavy meson remains as 
a spectator; for the 5c-meson decays there are two possiblities: h decays, while the c-quark 
remains as a spectator, and c decays, while the 6-quark remains as a spectator. In principle, 
in each spectator components there are two 'further' components the semileptonic one and 
the non-leptonic one: 



for 6-decay; 



Y[h^c)= E ^b^dv + E ^b-*cqq (6) 
l=e,ii,T q=u,d,s,c 



T(c^s)=Y: r,_r. + E ^c^sqq (7) 
l=e,fi q=u,d,s 



for c-decay. As for the concerned Be meson, being a double heavy meson, its two components 
6-quark and c-quark, each plays the decay role and the spectator role once in tern, so both 
of eqs.(^,|^ as the spectator components make contributions to Be decay. 

The semileptonic and non-leptonic decay rates of b— quark up-to the order 1/ml are 
evaluated by many authors Since in our numerical computations we need to use 

their formulas, so we quote them from the references into Appendix A. For c ^ s, the 
formulation is similar and even simpler, we also include the useful formulas in Appendix A. 



B. The non-spectator components in D and B meson decays 

The non-spectator contributions are crucially important to the D inclusive decays. For 
instance, the PI contribution may explain the data why r^ji ~ 2t£,o, but Tb± ~ t^o. More- 
over, the penguin contributions in D{Ds) and B{Bs) decays are negligible as aforementioned 
and the bound state effects emerge. Whereas, all the non-factorization effects cannot be re- 
liably well-determined yet. 
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With straightforward calculations, the precice operators for the non-spectator contribu- 
tions may be obtained. 

(a) For the D[Ds) decays: 

nif. ^ iV 

x[(l + - (1 + 2x^)B2] + 2c\[{l + ^)ei - (1 + 2a;+)e2]} 

2 2 

-ToVnsveJ^\Vcs?\Vus?Jl ' 4x+|(^ + 2CiC2 + iVc^) 

x[(l - x+)B, - (1 + 2x+)S2] + 2c?[(l - x+)ei - (1 + 2x+)e2]}, 

2 

r^^(D+) = ror/n,pec^|Kdr(|Kar(l - a;-)' + \VJ^) ■ [(c? + cl)iB, + 6ei) + 6ciC2i?i], 



ml 



2 2 
777- n iTT- lOiTT- lof/C' 



^^^(Z) + ) = -roVnspec^lKsl'Kdl'li^ + 2c,C^ + Ax2)(i?l - B2) + 20^(61 - 62)} 



-AT 



-row^iK.riK.r(i - ^+)'{[| + 2ciC2 + ivc^] 

X[(l + ^)i?l - (1 + 2x+)52] + 2C^[(1 + ^)6i - (1 + 2x+)62]}, 

r^^(D+) = ro77„,pec^|K.r(|V;.|'(l - x_)2 + \VX) ■ [(c? + c^)(5i + 6ei) + 6ciC25i], 

rW--.) = ^fc^^|K/(i--^)^ (8) 

oTT 771 n 



where 



^0 = T7^'Wec = 167r 

777 



1927^3''"-^- ^3 

2 



2^+ — ~Y']P+ — Pc + Pq: 

P+ 



— 2 

771: 



X-^ ^■,P-^Pc-Pq- (9) 



In the equations, the hadronic parameters are defined as follows: 



2mD, 



1 



2mn, '^^"^ ^ 8 

^ < d,\t:^id, > . te,, , (10) 



8 



where 



= n.T^Lqq^uT^Lc , (11) 



with T" = ^ and A" being the Gell-Mann Matrices, 
(b) For the B{Bs) decays: 

.2 

'N 



r^^(S^) = -ro77„,,ec|K<ir(l - z^r{{^ + 2ciC2 + Nclm + ^)B, - (1 + 2z+)B,] 



+2cl[[{l + f)ei-{l + 2z+)e2 



I f 

-ror^n^peclKdl Vl - ^^+{^N ^ + ^""^^f^^ ~ ~ + ^'"^^^'^ 

+2c2[[(l-^+)6i-(l + 2^+)62]]}, 

r^^(S-) = ro77„,pec^(l - -2-)'[(c? + C^)(Si + 6ei) + 6C1C2S1], 

r^^(S°) = -ro77„,,ec|K.r(l - z^f{{^ + 2ciC2 + Ncl)[{l + ^)B, - (1 + 2z^)B,] 
+2c?[[(l + ^)ei-(l + 2z+)e2 

-ror?n.pec|Ks| Vl - + 2CiC2 + 7Vc^)[(l - ^+)Si - (1 + 2z+)B2] 

+2c?[[(l-^+)6i-(l + 2;^+)e2]]}. (12) 



where 



^ — TT^TTT r c6 ) Vnspec — -l-OTr e— 

1927r'^ m? 



Analogous to D meson, the parameters Bi,B2, ei and 62 are defined: 

^'^"^ R in-? IR > -IkHl^n 

Bq 



£r < Bm > . (14) 
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where 



0% = h-i^Lqq-i^Lh , 



(15) 



C. The non-spectator components in B^. decays 

As pointed above, the spectator contribution to the Be hfetime is a sum of that from b 
and c individual decays as pointed above: 

Yspectator -psperfator _|_ Y^P^ctO'to^ (16) 

and r^P*^*^*"*"^ and Yf^^*"'*°^ are the same as they are in B and D decays and given in eqs.(|^J^). 
Now let us deal with the non-spectator contributions which are different from that in B and 
D decays. 

To estimate the non-spectator components in the B^ decays, let us write the relevant 
effective Lagrangian precisely here: 

L^f'f=\fi = = --^V,sV:,{c,{fi){s^,Lc){uYLd) + C2{fi){ui,Lc){sYLd)} + h.c, (17) 
and 

AG 



Lfff\fi = m,) = -_^{v;4K,(ci(/i)0r + C2(/i)0^) + K:(ci(/i)OJ + C2(/i)0- 



+ 



^ h,LucYLb + V:,Y.Ci04} + h.c., (18) 

l=e,T,ii i=3 



where the operators are 



01 = sj^Lcc^Lb, 
01: = d-f^Luc-f^Lb, 

02 = Si'j^LcjCj'j^Lbi, 
O2 = di'j^LujCj'j^Lbi, 
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Os = s'jf.LbcYLc, 
O4 = Si'y^LbjCjYLci, 

Oe = Si'y^LbjCjYRci 



(19) 



and Ci{i = 1, 2, ■ ■ ■), denoting the Wilson coefficients due to QCD corrections, will take the 
values as those in ref. |Q. Here we consider the non-spectator components in Be decays by 
two steps. The ffist step is to compute the relevant operators upto the order OII/itlq) and 
then to evaluate the contributions precisely. 

1. Pauli interference (PI) operators 



The Pauli interference (PI) operators P^gg and TpJ^g^^^ which correspond to the non- 
leptonic decay induced by the tree part and penguin respectively are given by: 

2Gl 



tree 



penguin 



where 



TT 



\2 2 



{2ciC2 ■ b%Ldc>-f>'W + (c? + cl) 



■b'-f^Lc^c>YLb'} 
2G^ 



\Vcb\'^\Vcs\'^{ 



1 - Z.fpl ■ { (2ciC3 + 2C2C4 + 2C3C4) ■ U^^Ldc>^^W 



+ {cl + cl + 2ciC4 + 2C2C3) ■ U-i^Lc>c'-i^LU] 



G 



+^W\V,s\\l - z.f • {(1 - z.)v'_g^^ + 2(1 + 2z.)jf_f\ 



■{2C5C6 ■ b%Wc^^,Rc' + (4 + cl) ■ b%m^,Rci} 

J TTigP • |[C2C6 + C3C6 + C1C5 + C4C5J 



IKriKsi 



■[U-i^'Ldc>-ic,-if,W + b'-ff,-faRc"c^-f''Lb>] + [C2C5 + C3C5 + ciCg + c^cq] 

■ [¥rLc^(?lal,Lb' + Fj^j^Rc^c^rm } , (20) 



— 2 



2 ' 



P 



P- 



Pb-Pe 



(21) 
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2. Weak annihilation (WA ) operators 



The weak annihilation operators are T^^, ^^enguin ^"^^ f ^^(5^ — tu^.) which correspond 
to the non-leptonic decay induced by the tree part, penguin and the pure leptonic (PL) decay 
respectively.0 



tree 



penguin 



■{{Ncl + 2ciC2) ■ U^^Ld(9^,LV + cl ■ 



(22) 



37r 



W+} 



■{{Ncl + 2ArciC4 + 2C3C4 + 2ciC3 + 2C2C4) ■ V-i^Ldc^-i^LV 



Vcb?\Vcs\^{l - z+)2p^ • [{Ncl + 2C5C6) ■ URddW + ■ VRddLU 



TC 



P(l - z+fm 



,P1 ■ {[Ncq{Ci + C4) + C5(Ci + C4) 



+C6(C2 + C3)] ■ [b'Rdc'-f^W + 6'7^Lc*e'L6^] 
+(C2C5 + C3C5) ■ [URdc'^^LU + 6*7^Lc^c'L6^]}, 



26-2 



cb\ 



37r 



{{l + "-^)pl9'^-{l + 2z^)p%p\} 



(23) 
(24) 
(25) 



where the parameters p+, 2+ and 2,- are defined by 



P+=Pb+ Pe 



pi 

2 

P\ 



ml 



(26) 



^Due to hilicity suppression, the decays Be — > l{e,fj,) + are ignorable for the low order estimate 
of the lifetime, thus we do so here. 
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3. The contributions from the non-spectator WA and PI to the lifetime for Be meson 



With the optical theorem, substituting all the above operators f ^^^^ f into the relevant 
matrix element, we may estimate the non-spectator contributions to the lifetime of meson: 

1 



2Mf 



-{Bc\T\Bc^ 



(27) 



where f denotes the relevant operators for PI and WA given in the above subsections. 

According to eq.(P7| to evaluate the lifetime, finally some hadronic matrix elements need 
to be determined, whereas, having nonpertubative nature, they cannot be determined by 
well-established theories so far. Let us discuss their deturmination here for a while. 

First of all, some parameters, such as -B — 1, B2, Bi, B2,ei, e2,ei and €2, appear in the 
corresponding estimates for B and D decays too. Precisely for Be decays, they are 

1 



2Mb, 
1 

2M^ 
1 

2M^ 
1 

2M^ 
1 

2il^ 
1 

2M^ 
1 

2M^ 
1 

2Mr 



{Bc\Ol._ 


-a\Bc) 


{Bc\Ol 


-p\Bc) 




-a\Bc) 


-{Bc\Tl 


-p\Bc) 


{Bc\d'y_ 


-a\Bc) 


{Bc\d^s- 


-p\Bc) 


{Bc\f^_ 


-a\Bc) 


-{Bc\fl 


^p\Bc) 



PbcMb. 



PbcMb. 



PbcM^ 



Be 



PbcMb. 



Bi, 

B2, 

-ei, 

-62, 



8 ^' 
PbcMb. 5 

7^ ^2, 



where the relevant four-quark operators are 

Oy_A = hf^LccYLb, 
Og_p = bLccRb, 
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-p 


= bLT''ccRT%, 




-A 


= Irf^RccYRb, 




-P 


= bRccLb, 




-A 


= hf^RT''cc^^RT%, 




-P 


= bRT^ccLT^b . 



(28) 

There are eight extra matrix elements corresponding to the newly emerged operators in the 
Be case. The 'new' matrix elements relate to the above parameters or new ones (es, e^, €5, ee) 
as follows: 

^ {Bc\bLccLb\Bc) = zlhM^B^, 



2M, 



Be 



^ {Bc\bRccRb\Bc) = — %^S2, 



Be 



2M, 



2M 

^ {Bc\h>.LccYRb\Bc) = 



Be 



^{Bc\h,Rcci^Lb\Bc) ^ 

{Bc\bLT''ccLT%\Bc) = •' ^^ 63, 



2M, 



Be 



2M, 



^ {Bc\bRT^ccRT%\Bc) = -J^^^e^, 



Be 



2M 



^ {Bc\b-ff,LT''ccYRT%\Bc) = — %^e5, 



Be 



1 _ fl Ji/T 

/D„II,., TDH^CL -^LL T rria-l\ TZ) „\ J Bc ' 



{Bclbj^RT^ccj^LT^lBc) = -^^^f^e^. (29) 



2Mb, 

Generally speaking, we may assume that 

Bi(2) = Bi{2)] ei(2) = ei(2), (30) 

with symmetry consideration. As for the parameters €3,4 and e^^, we would conjecture that 
€3,4 ~ €2 and £5,6 — ei instead of precise computation^. 



^The conjecture should be tested and proved later on. It should be considered as a working 
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In the earlier literatures, usually Bi ^ B2 ^ 1 and ei ~ —0.15 from the lattice calcula- 
tions and 62 = 0. According to our numerical computations and trials to fit the data about 
the lifetimes of the heavy mesons D^, D^, Dg, B^, B^ and Bs and their semileptonic decay 
branching ratios as well, we find that to adjust the values of them and the pole masses of b 
and c quarks, when the parameter £2 7^ etc, a better fit is obtained. 

Now for the nonspectator component PI, we have 

Kie = ^fLMsM,m.\\l - z.)V- ■ {[2C1C2 + ^(c? + cl)]B, 

+2(c? + c^)ei} (31) 

^Z9u^n = ^PbcMbMKsWI ' Z-)Y_ ■ {[2C2C, + 2c,C, + 2c,C, + (32) 

1 

iV 



j^i4 + 4 + 2C2C3 + 2ciC4)]Si + 2(4 + cl + 2C2C3 + 2ciC4)ei} 

5/LMBjK.riK/(l - z.r . {[2c5Ce + Ucl + cl)][^^p^B2 



Gj 

^ ^ {mlBi + mlBi - 4mbmcB2 + 2mbmcBi)] + 2(c^ + c^) [^-^^p?.e2 



6 

"^5 ^~ (^ft^i + "^c^i ~ 2mbmc{e3 + £4) + mbrndes + ee))]} 

--^fBcMBjVcb\'^\Vcs\'^{l - Z-frflc ■ |[ciC5 + C2CQ + C3C6 + C4C5 
OTT 

+-^(ciC6 + C2C5 + C4C6 + C3C5)][2mc-Bi + mb{-4:B2 + 2Bi)] 

+2(ciC6 + C2C5 + C4C6 + C3C5)[2mcei - 2mfe(e3 + £4) + mb{e5 + ee)]} , (33) 



and for WA, we have 



- -^\Vc,\WfLMB.{l - z^f ■ {[Ncl + 2c,C2 + |] (34) 
x[(l + ^)Ml^B, - (1 + 2z+){mlB2 + 771^^2 + 2mbm,B2)] 
+2cl[{l + ^)Mlei - (1 + 2z+)(m^e2 + mlh + nibm^ie^ + 64))]}, 

r^^(i?e - TV) = -^WKs?fLMsXl - z^f ■ {(1 + \)MlB, (35) 



assumption. With the assumption eq.(30) in addition, all the assuptions here will cause an essential 
uncertainty for the estimates. 
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-(1 + 2Zr){mlB2 + mlB2 + 2mhm^B2)}, 

Ken^mn = " ^ I H IVIM^, ( 1 - Z^)' ' {[(^^^ + + C,){cs + Nc,)] (36) 

x[(l + y)M2^5i - (1 + 2z+)imlB2 + mlB2 + 2m,m,B2)] 

+2(2c2 + C3)C3 ■ [(1 + yHei - (1 + 2z+){mlt2 + mlh + rrn,m^{e^ + 64))]} 

+^|K=6riK.| VlcMi(l - z+f ■ {[§ + 2C5C6 + N^^B2 + 20^62} 

-^|V;,nK.|VleMB.m,(l - Z+f ■ {[(^^ + Ci + C4)(C5 + N C,)] 

x[2mbB2 + 2mc52] + 2(c2 + 03)05 ■ [mb(e3 + £4) + 2mfce2]}- 



D. The effective mass of the decaying heavy quark 

The masses of the acting heavy quarks in a decay must be treated carefully although the 
bound-state effects make the problem complicated and obscure. It is commonly accepted 
that if the charm quark appears as a decay product, the mass should be its running one at 
the energy scale of the decaying quark or the meson, whereas, if it appears as the 'parent (s)' 
of the decay, the quark (antiquark) is not "free", but in a bound state, thus the pole mass 
should be taken and the bound-state effects on the mass must be taken into account too. 
Especially in spectator mechanism the decay possibility of the heavy quark is very sensitive 
to the value of its 'adopted' mass, hence what a value of the quark mass adopted in the 
estimate must pay special care. Narison [IIH] used the QCD sum rules to estimate the mass 



difference — My^-^ where M^^"^ is the short-distance perturbative pole mass and M^^ 
is the long-distance QCD-related effective mass up-to two-loops. Whereas the authors of 
||T7| attributed such effects into a factor which is multiplied to the decay width of the "free" 
quark. 

Here instead of deriving the modification factor with a relatively large uncertainty, we 
treat the problem phenomenologically i.e. by introducing a parametrization 

M^^^ = Ml°^' - A, (37) 
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where A manifests the bound-state effects, and it will be fixed phenomenologically. Note 
here that for each heavy meson there are three quantities: lifetime (total width), inclu- 
sive semileptonic branching ratio and pure leptonic branching ratio which may be used for 
phenomenological analyses, so the estimates here are still well-determined even when we 
introduce the parameter A here. 

In the next section, we will discuss A and other related parameters more precisely. 

With all the formulae derived above and the hadronic matrix elements, we can make 
numerical evaluation of the lifetime of Be straightforwardly. 



III. NUMERICAL RESULTS 



Since we carry out the estimate of the lifetime of Be with a 'gloable' comparison to all of 
the heavy and double heavy mesons, so the determination of all of the parameters by fitting 
the existence experimental data is 'over-determined' for our goal and has certain level tests. 
Therefore we evaluate the lifetimes, the semileptonic branching ratios and the pure leptonic 
branching ratios for all the mesons D^, D^, D^, B^^ 5°, Bg and B^. in this section in turn and 
present the numerical results in this section. 



A. For the heavy mesons D and B 

To evaluate the lifetimes of D°,D'^,Ds,B°,B^, 5, mesons and their branching ratios of 
the semileptonic decays, we use the formulae given in sections 2.1, 2.2 and the appendix. 
The parameter values are taken as follows. \Vcs\ = 0.974, \Vud\ = 0.975, as{mc) = 0.29, 
ci{mwc) = 1.30,C2(mc) = -0.57 [g, Bl = B2 = l,ei ~ -0.05, = §, the decay 
constants of D mesons = 160 MeV, fn^ = 190 MeV. In the evaluation of the Pauli 
interference PI contribution to D decay width, we take the = (pc ~'PqY value as 0.5 M|) 



as done in ref. [p^] . 
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We take mf'*" = 5.02 GeV, mf'^ = 1.88 GeV By eq.®, we have the running 



mass of the charm quark at various energy scales as 

mdmc) = 1.67 GeV, mdmb) = 1.41 GeV, m^{mB,) = 1.37 GeV. 

By fitting data [Q, we should have the quark masses as = 125 MeV, mj;-^-^ = 1.65 GeV 
respectively. Then we obtain the D meson lifetimes: t{D^) = 0.419 ps, t{D~) = 1.06 ps, 
t{D~) = 0.446 ps, and the branch ratio of the semileptonic decay of meson Bsl{D^) = 
6.9%. Comparing to the experimental data: r(D°) = 0.415±0.04 ps; t{D^) = 1.057±0.015 
ps; t{Ds) = 0.467 ± 0.017 ps and Bsl{D^) = 6.75 ± 0.29%, one can see the fit is quite well. 

For the estimate of B meson lifetimes, we take the mass of the charm quark at final 
states rric to be running mass, namely it is different from the pole mass of charm quark 
in D-decays but the running mass at the energy scale mh i.e. mc{mii) = 1.41 GeV. When 
calculating PI contribution to the B~ decay width, we take the value of = {pb — PuY 
approximately to be 0.8 [^]. For the other parameters in the numerical computations 
the values are adopted: \Vcb\ = 0.04, asijnb) = 0.20, Ci(mfe) = 1.150, C2{mb) = —0.313 0, 
a = 1.06, (3 = 1.32 Bl = B2 = 1 and d = -0.14, 62 = -0.08 [||]. The decay constants: 
Jb = 200 MeV and Jb^ = 220 MeV. Furthermore, taking the b quark pole mass m^"'^ = 5.02 
GeV and ml^^ = 4.89 ~ 4.91 GeV, we obtain the results that 

r(fi°) = 1.54ps, t{B-) = 1.74ps, t{B;) = 1.56ps, B,i{B°) = 11.2%, if m^^^ = 4.89GeV 

r(5°) = 1.52ps, t{B~) = 1.71ps, t{B;) = 1.54ps, 5,^(5°) = 11.2%, if ml^^ = 4.90GeV 

t(S°) = 1.50ps, t{B-) = 1.68ps, t{B;) = 1.51ps, B,i{B^) = 11.2%, if ml^^ = 4.91GeV 

where Bgi indicates the branching ratio of the semileptonic decay. Comparing with the 
experimental data r(5°) = 1.56 ± 0.04 ps, r(5±) = 1.65 ± 0.04 ps, r(5,) = 1.54 ± 0.07 ps 
and Bsl{,B^) = 10.5 ± 0.008% we can see the fit is quite good. Accouding to the definition 
of A, we have = mP°^^ - miff = 0.23 GeV and A5 = - miff = 0.11 ~ 0.13 GeV, 
that is understandable. 
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With all the parameters obtained by fitting the lifetimes of B^, B^, Bs, D^, D^, Dg and 
the branching ratios of the semileptonic decays of B and D mesons, we are proceeding to 
evaluate the lifetime of i?c— meson and its semileptonic decay rate. 

B. For the double heavy meson Be 

The spectator componet contribution to the Be lifetime is a sum of the individual b 
and c quark decays, while leaving the other one as a spectator. When evaluating this 
contribution, mi, is its pole value at = ml, and m^ also the pole value mdmc). Whereas 
for the non-spectator contributions, i.e. the WA and PI pieces, the corresponding energy 
scale for the running charm-quark mass in the final state is taken as Mb^- Now let us 
take the relevant parameters for B^. as follows: Mb^ = 6.25 GeV, M^^ = 6.33 GeV, Bi = 
B2 = l,ei = —0.14,62 = —0.08. For the decay constant, we adopt Eichten and Quigg's one 
fsc = 500 MeV and the lattice one Jb^ = 440 MeV respectively. Furthermore in 
the calculation of the PI contribution, the quantity = [pf, — pc)"^ ^ 2ml + 2m^ — M^^ is 
taken approximately. With the parameters described above, we obtain the numerical results 
and tabulate them in Table 1. 



Table 1: The results for B.. meson 







'ppen. 






YWA 


YPI 




BsL 


440MeV 


0.362 (ps) 


3.4% 


22.8% 


70.9% 


13.4% 


-7.1% 


0.078 ps-i 


8.7% 


SOOMeV 


0.357 (ps) 


4.3% 


22.4% 


69.7% 


16.9% 


-9.0% 


0.100 ps-i 


8.4% 



In the table /^^ denotes the decay constant; tb^'- the lifetime of Be] F^^"': the contribution 
from the interference between the penguin and 'tree' terms; F(rz/): the width of the pure 
leptonic decay (r chennal only but almost equal to the total) and the Bsl- the branching 
ratio of the semileptonic decay of the meson B(.. 

Since both the 'parents' h and c quarks reside in a bound state {B^. meson), the problem 
how to choose the value of the masses mi, and m-c emerge as in the cases of the heavy mesons 
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D and B etc, but when taking all the parameters fixed by fitting data as the above, we 
obtain the results presented in table 1. 

Let us discuss the bound-state effects on the b and c quark-masses in Be meson more 
precisely. Because Be includes two heavy quarks i.e. is a double heavy meson, the bound- 
state effects might be greater than in the heavy mesons B,D. We think that the values 
m^-^-^ and ml^^ might be smaller than m^-^-'^ = 1.65 GeV and ml^^ = 4.9 GeV that we 
obtained in B and D decays. Phenomenologically, if in Be meson, rrfj^ [Be) = 1.55 GeV, 
ml^^ (Be) = 4.85 GeV, we obtain T{Be) ~ 0.47 ps, which occationally is closer to the center 
value of the Be lifetime measured recently |l|. In this case, Ac = 0.33 GeV and A;, = 0.17 
GeV. Because the rates of direct h and c decays, which dominate the lifetime of Be meson, 
are proportional to (Mq-^)^, the results are so sensitive to the effective masses. 



IV. CONCLUSION AND DISCUSSION 



In this work, we estimate the lifetime of Be in the 'unique' theoretical framework where 
the nonspectator effects are taken into account properly and the necessay parameters are 
determined by fitting the data of the heavy mesons 5°, B^, Bs, D^, D^, Ds on the lifetimes 
and the branching ratios of the inclusive semileptonic decays as input. 

Not only the uncertainties in the estimate are discussed but also the physical parameters 
appearing in the estimation are fixed in reasonable regions. 

Even though not all of the parameters are fixed by fitting the avialable data for the heavy 
mesons B^, B^, Bg, D^, D^, Dg, in order to carry out the estimate we make some reasonable 
assumptions or conjections. In fact, in terms of the lattice calculation, QCD sum rules and 
other approaches, the parameters, such as Bi and B2 (the factors in the hadronic matrix 
elements, the values manifest the deviation to the vacuum saturation), may be determined 
at certain accurate level. The deviation from model-dependent calculation is small, and 
can be neglected in practice. Some of the parameters, such as ei and €2, being realized to 
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relate to the non-factorization effects |33], may be calculated at the D-meson scale and the 



B-meson scale respectively in terms of the QCD sum rules P, p8| , p9| . It is known that the 
numerical values obtained by the QCD sum rules may have errors about 10 ~ 15%, but they 
still can be used as inputs to the phenomenological calculations without causing too large 
errors. In this work we also carefully consider the quark masses and take into account of the 
bound-state effects. The consistency of our numerical results of B and D mesons confirms 
the validity of the parameter regions. 

The earlier estimations on the lifetimes of B and D mesons and the semileptonic decay 
rates obviously deviate from the data. Luke, Savage and Wise pointed out that in decay 
c Xez/g, the contribution of a'j. order is of the same magnitude as that of 0{as) and 
this higher order correction suppresses the semileptonic decay rate of D-meson. Taking into 
account this fact, we obtain numerical results of lifetimes of D-mesons and their semileptonic 
decay rate, and find that they are satisfactorily consistent with data. Whereas, for B-meson 
decays, the order correction, as well as the 0{as) correction, are smaller. With these 
corrections concerned, the results for B-mesons are also consistent with data within the 
experiment tolerance region. All these imply that the parameters taken as the above are 
reasonable. 

When evaluating the Be lifetime and its inclusive semileptonic decay rate, some new 
aspects must be taken into account. First there are several new operators in the effective 
Lagrangian playing roles. Their appearance is due to non-negligible m^, whereas in B and 
D cases, the light quark mass nig is ignored with quite high accuracy. Correspondingly, 
several new hadronic matrix elements are induced by these operators. Some of them are 
also proportional to Bi and B2, which exist in the expressions for B and D meson decays, 
as long as the factorization theorem and the vacuum saturation hold. But in the non- 
factorization contributions, new parameters appear and 63 ~ eg is assumed. In this work, 
we have taken a naive symmetry consideration to let €3,4 ~ €2 and €5^ ^ ei. 

As discussed in the introduction, in the B^ case, the interference between the penguin 
and tree terms is not negligible. Namely, the penguin contribution to Be lifetime is much 
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more important than to B and D decays. Our results confirm this allegation and we have 
found the contribution from the interference can be as large as 3 ~ 4% of the total width. 
This fraction is measurable in accurate measurements. Since direct measurement of penguin 
diagrams would be interesting, these sizable value can be encouraging for future experiments. 

The lifetime of i?c-nieson is estimated as 0.37 ~ 0.367 ps for /^^ ~ 440 ~ 500 MeV, which 
are smaller than the central value of the measurement tb^ = 0.46to;i6(stat)±0.063(syst) ps 
Our earlier estimation of tb^ @ was 0.4 ps. Anisimov et al. estimated tb^ in the 
light-front constituent quark model and obtained tb^ = 0.59 ± 0.06 ps, which is larger than 
the measured value. In our estimation, we use the values of Mb^ and M^^ as 6.25 GeV and 
6.33 GeV 0, whereas the measurement is Mb, ~ 6.40 ± 0.39(stat)±0.13(syst) GeV. When 
the bound-state effects on b and c masses are reasonably taken into account, we can have 
tBc ~ 0.47 ps, which is very close to the present experimental center value of the i?c-nieson 
lifetime. As noted, the change of /b, itself only does not influence the result much, e.g. 
as J'bc changes from 440 MeV to 500 MeV, tb, varies 1% only. Whereas, a large error in 
mass of the decay parents would result in an estimate error of about 4%, therefore, the more 
accurate lifetime and mass of meson will test the framework adopted here the deeper. 

The more accurate experimental measurements will shed fresh lights on the framework 
adopted here, which involving the effective heavy flavour theory and the duality between 
quark states and hadronic states 

i,j k 

etc. It will help to clarify many uncertainties. 
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Appendix A 

The semileptonic and non-leptonic decay rates of b quark through order 1 / are given 
as following P,p5|. 



mt 



/l(Xe,0,0) 



Tr^LiH,) = r[,') ■N-{{cl + cl + ^) ■ [(a/o(xe, 0, 0) + (3Io{x,, x,, 0)){H,\bb\H,) 



mi 



N 

'MXc,0,0)+h{Xc,Xc,0))] 



mf^ N J 



(38) 



where 



in = , _ _ „ I Vcb\ 



° 1927r3 



(39) 



and the following notation has been used: Iq, Ii and I2 are phase-space factors, namely 



/o(x,0,0) = (1 - a;2)(l - 8x + x^) - 12xMoga;, 
Ji(x,0,0) = ^(2-x^)/o(a;,0,0), 
/2(x,0,0) = (l-x)3, 

Ux, X, 0) = v(l - Ux - 2x^ - 12x^) + 24a;2(l - x^) log 

\ — V 

/i(x,x,0) = i(2 - x-^)/o(x,x, 0), 
2 ax 

hix, X, 0) =v{l + ^ + 3x^) - 3s(l - 2x^) log ^ ^ ^ 



1-v' 



Xc = {mjmbf, t; = Vl - 4x , 



(40) 



with Io^i^2{x,x, 0) describing the b — > ccs transitions. 

And for ri{xc, xi, 0), which is the QCD radiative correction to the semileptonic decay rate. 
Its general analytic expression is given in [^. The special case ri{x, 0, 0) is given in |^ and 



it can be approximated numerically by p2,25 
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r/(x,0,0) 



2a. 



>^-y)(l-v^)^ + | 



(41) 



For the decay b cru, according to [0] we roughly have 

T{b CTu) ~ 0.25r(6 ^ ceu). 
The expressions are simpler for c — > s: 

Tsl{H,) = ri;"^ ■ r]{xs,xi,0)[lo{xs,0,0){Hc\cc\H,) 

2CiC2, 



(42) 



2(/i|) 



/i(a;.,0,0) 



r;vL(i^c) = ■ AT . {(c? + + ^) X [a/o(x„0,0)(i/,|cc|i/,) 



'-^^hix^, 0, 0)] - 8i%^ . 0, 0)}. 

mi mi i\ > 



(43) 



(44) 



where 







192^3 1^-1 ' ^« 



— 2 



(45) 



and for the correction ri{xsi 0) in the c-decay case, we adopt a numerical expression from 
lEIl. It reads 



risL 



1 -2.08( 



n 



22. 7( 



TT 



(46) 



For the dimension-three operator QQ, the expectation value can be expressed at follows: 



{Hq\QQ\Hq) = 1 



2ml 



+ 



+ 0{l/ml) 



(47) 



where ((pg)^) = {HQ\Q(iDfQ\HQ) denotes the average kinetic energy of the quark Q 
moving inside the hadron and {ij^'g)hq = {HqlQ^a ■ GQ\Hq). 

Based on refs. the kinetic terms take the values respectively as follows: 



Be n n/1 ((Pc)^)Bc 



0.016 , 



0.21 



0.04 



0.4 



''b '"'b 

For the chromomagnetic operator one finds {fJ'G)pQ — ^mq^My^ — Mp^), where Pq and Vq 
denote the pseudoscalar and vector mesons, respectively. 
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(48) 
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